Imaging density disturbances in water with 41.3 attosecond time resolution 
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We show that the momentum flexibility of inelastic x-ray scattering may be exploited to invert 
its loss function, alowing real time imaging of density disturbances in a medium. We show the 
disturbance arising from a point source in liquid water, with a resolution of 41.3 attoseconds (4.13 x 
10~ 17 sec) and 1.27 A (1.27 x 10 -8 cm). This result is used to determine the structure of the electron 
cloud around a photoexcited chromophore in solution, as well as the wake generated in water by 
a 9 MeV gold ion. We draw an analogy with pump-probe techniques and suggest that energy-loss 
scattering may be applied more generally to the study of attosecond phenomena. 



Brisk progress has been made recently in the gener- 
ation and detection of ultrashort, attosecond (1 as = 
10 -18 sec) laser pulses with high harmonic generation 
techniques H E B II IE II 011 This has heralded 
an age of attophysics in which many-electron dynamics 
will be probed in real time. Much of what is already 
known about electron dynamics has been derived from 
energy-loss techniques such as inelastic x-ray, electron, 
or neutron scattering, which have been fruitfully applied 
to the study of, for example, plasma oscillations, exci- 
ton dynamics, and spin waves. Such experiments are 
normally examined in the frequency / momentum repre- 
sentation, where the spectra may be easily compared to 
theoretically calculable response functions for the elec- 
tron density or local magnetic moment. However if such 
measurements truly reflect dynamics, a temporal repre- 
sentation may also be illuminating. 

In this letter we demonstrate a method for inverting 
energy-loss measurements into time and space, permit- 
ting explicit imaging of electron dynamics in a medium. 
This allows the spatial extent of excited states to be de- 
termined, i.e. the mean free path or dephasing distance, 
and a parallel with pump-probe techniques to be drawn. 
Inverting requires reliable sampling of a loss function 
along both energy and momentum axes with sufficient 
range to resolve the features of interest and sufficient res- 
olution to provide an adequate field of view. Therefore, 
because of its kinematic flexibility, we consider the case 
of inelastic x-ray scattering (IXS). 

In IXS a photon with well-defined initial momentum 
and energy, (ki,cjj), is impinged on a specimin which 
scatters it to a final (k/ , Uf) . The spectral density of scat- 
tered photons is proportional to the dynamic structure 
factor of the material, S(k, to), where u) = u>i — LOf and 
k = k; — k f are the transferred energy and momentum 0. 
S (k, lj) was posed originally by van Hove[I3 as a mea- 
sure of the dynamical properties of an interacting elec- 
tron system, and by the fluctuation-dissipation theorem 
is related to the imaginary part of a response function |ll| 

Im [ X (k, w)] = -7T [S(k, u>) - S(k, -w)] (1) 
which can thereby be experimentally determined. x(k, u) 



is known as the density Green's function, or density 
propagator, and describes the way disturbances in the 
electron density propagate in the system. %(k, u>) is 
the space-time fourier transform of x( x : t) = —i < 
0|[5n(x, t), Sh(0, 0)] 1 > 6(t), where 5n is the density fluc- 
tuation operator, and in real space represents the distur- 
bance produced by a delta function source at the origin 
at t = 0. Energy loss scattering, in this sense, is rather 
like a pump-probe experiment; the system is perturbed 
at a reference point in time, and its subsequent evolution 
observed. To make the analogy more concrete we wish 
to carry out such an inversion on a real IXS data set. 

Unfortunately, inverting requires knowledge of the full 
x(k, u>), but the experiment provides only its imaginary 
part. In other words in IXS there is a phase problem, 
like that in x-ray crystallography, that must be over- 
come before results can be visualized explicitly. For- 
tunately, the phase can be retrieved by exploiting the 
causal properties of x(k, u>), i.e. the fact that it satisfies 
the Kramers-Kronig (KK) relation, the second of which 
may be written [l2| 
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Eq. (2) ensures that x( x , t) =0 for t < and is equivalent 
to choosing a retarded causality convention for the prop- 
agator. From (2) Re[x(k, u)] may be determined, which 
in principle allows reconstruction of y(x. £)[l3l ITij. 

We used this procedure to image the density distur- 
bance generated by a point perturbation in liquid water. 
The elementary quantum of density oscillation in water 
is the plasmon, which arises from collective vibration of 
the 2s and 2p shells of Oxygen whose density determines 
the normal frequency ui p = 22 eV jl6t| . The plasmon is 
also the fundamental source of electronic screening in 
water and determines, for example, its optical refractive 
index[l7| . We used IXS to sample the plasmon momen- 
tum and energy dependence over the ranges 1 eV < uj < 
100 eV and 0.476 < k < 4.951" 1 (Figure 1). Mea- 
surements were carried out on station C-l at the Cor- 
nell High Energy Synchrotron Source (CHESS), where a 
nested Si(lll)/(400) monochromator with sagittal focus- 
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ing was employed [18j, and at CMC-CAT at the Advanced 
Photon Source (APS), where a single Si(333) channel cut 
was used. Both experiments used a diced, backscattering 
Ge(444) analyzer to achieve overall energy and momen- 
tum resolutions of Aoj = 0.30 eV and Ak = 0.146A" 1 . 
This configuration provided effective spacetime resolu- 
tions of At = 41.3 as and Aa; = 1.27A, with fields of 
view of 13.8 fs and 54.3 A. 

In addition to inelastic scattering the spectra contained 
zero loss lines due to the static liquid structure factor. 
This was subtracted and the spectra extrapolated to at 
oj = 0, where 7ro[x(k, oj)] must vanish by causality, and 
at k = 0, where it must vanish by charge conservation. 
The F sum rule was applied with no = 0.20/A 3 to place 
the spectra on an absolute scale. 

While the proposed Kramers-Kronig inversion proce- 
dure is simple in principle, unexpected subtleties arise 
when applying it to a discrete and finite data set. First, 
a discrete Im[x(w)] implies that x(t) is periodic, a prop- 
erty incompatible with the causality constraint that x(i) 
vanish for all t < 0. Second, the KK relations are de- 
fined on an infinite oj interval, yet our scan range was 
a finite 100 eV. We resolve the former issue by analyti- 
cally continuing the spectra onto a continuous frcqcncy 
interval, by linear interpolation, resulting in a time axis 
which is formally infinite. This permits to vanish for 
all t < 0, but brings about an aliasing effect that makes 
the dynamics repeat with a period T = 2-k / doj — 13.8 fs, 
where doj is the spectrometer resolution (here the quan- 
tity 2ir/doj plays the role of the Nyquist critical frequency 
familiar from Fourier optics). 

The latter issue is resolved by extrapolating the spectra 
to oj — » oo with a Lorcntzian fit, which is appended to 
the measured spectra (Figure 2). This forces the time 
domain to be formally continuous, causing x(x, t) to be 
defined even on times much shorter than the effective 
resolution of 27r/100 eV = 41.3 as. It is fair to examine 
the function on such short time scales, though one must 
remember that its behavior is sensitive to the form of the 
extrapolation (Figure 2). 

Figure 3 displays time frames of x( x , i). Because the 
medium is isotropic the disturbance is spherically sym- 
metric, and each image is a section of this sphere with 
X(x, t) plotted in units of A~ 6 on the vertical axis. At 
negative times x(x, t) vanishes and the system is "placid" 
(Fig. 3a). At t = the system is "struck" with a posi- 
tive perturbation, generating a negative recoil at the ori- 
gin surrounded by a positive build-up at |a;| = 1.3A as 
current flows away from the source point (3&). On an 
expanded scale a dip is also visible at |x| = 2. 5 A, fol- 
lowed by another peak at 3. 8 A (3h), forming a transient 
structure analogous to the Friedel oscillations that occur 
around point impurities in metals [l^.l20| . 

Because of the compensating charge from the ion cores 
the pattern (3/i) experiences a return force which causes 
it to reverse direction and change shape (3 c). An anhar- 



monic oscillation occurs, with a time scale of order the 
plasma frequency T — 2tt/oj p — 180 as (3d-g). Over- 
all, the disturbance is seen to be rather local, occurring 
within 5 A of the origin, a property related to the broad 
line shape in Figure 1. A plasmon in water can oscillate 
a few times but never develops into a true propagating 
mode. 

As time evolves the entropy grows and the disturbance 
damps as the system evolves toward thermal equilib- 
rium. The disturbance decays below our threshold of 
dctectability, determined by the experimental noise level 
shown in Fig. 3z, after an elapsed time of 350 as — less 
time than it takes light to travel 100 nm in vacuum. 

x(x, t) provides a means to visualize electron dynam- 
ics, but more significant is its relevance to a rather broad 
class of ultrafast processes in matter. Specifically, be- 
cause it describes the disturbance from a point source, 
by the principle of superposition x(x, t) can be used to 
determine the effects of extended sources. The charge 
nindfat) induced in a medium by a time-dependent 
source n ex t (Xl£) is determined in reciprocal space by the 
relationship 2 lj 



n iri <i(k, oj) = -^-x(k, oj)n ext (k, oj) (3) 

which in real space is a convolution over the dimensions of 
n ext (x, t) identical to the Green's function integral used 
to handle sources in ordinary differential equations |22|. 
We illustrate its use by modeling an oscillating dipole, 
and a gold ion travelling in water at 0.01 times the 
speed of light, c. As a model of the former we take 
n ex t(x, t) = [S(x + 0.5 A) — S(x — 0.5A)] cos(woi), i.e. two 
point charges with opposite sign separated by 1 A, os- 
cillating with a frequency ojq — 2eV/h. The resulting 
disturbance, n ext (yi,t), is shown in Figure 4a. Such a 
disturbance would occur, for example, around a pho- 
toexcited chromophore in solution, and should have influ- 
ence on short-range fluorescence resonance energy trans- 
fer (FRET) commonly used to measure distances in bio- 
logical systems 23J . Similarly, an ion with charge Z trav- 
eling at velocity v, represented as n ext (x, t) = Z5(x— vt), 
for Z = 79 and v = 0.01c produces the pattern in Figure 
4 b. This is an image of the density wake of a 9 MeV gold 
ion traveling in water and can be used, for example, to 
quantify the stopping power, dE/dx[24^. 

It is important to address the issue of why, despite 
the spatial resolution of dx = 1.21 A, individual atoms 
are not visible in Figs. 3—4. The reason is that a den- 
sity propagator is actually a function of two spatial vari- 
ables, x( x , x';t), with a Fourier transform of the more 
general form x(k, k';w). There being only one momen- 
tum transfer in scattering, we measure only the diago- 
nal (i.e. longitudinal) response, x(k, k;w) 11]. In other 
words the present analysis has assumed a translational 
invariance that does not strictly apply. The images in 
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FIG. 1: False colour plot of — 7?n[x(k, a>)] for water in units 
of as/ A 3 , plotted against transferred momentum (horizontal 
axis, in A -1 ) and energy (vertical axis, eV" 1 ). The broad 
feature is the valence plasmon, corresponding to collective 
oscillation of the 2s and 2p shells of oxygen. 



FIG. 2: Individual spectra at k = 0.583A -1 (above) and 4.95 
A -1 (below), showing the raw data (open circles), extrapo- 
lated tails (red lines) and the resulting real part (green) line. 



Figs. 3—4 must therefore be thought of as spatial aver- 
ages — not single events but the mean of many events 
at different locations in the specimin. In many cases an 
ensemble average provides adequate information; if not 
this limitation can be overcome by using standing wave 
techniques [l^ . 

The existence of nonlocality in electrodynamics, in 
both time and space, is well established |l2j, |2j|. I n this 
article we have shown that they may be combined to ex- 
plicitly image electron dynamics. Because x-rays may 
be tuned over extremely broad ranges of energy, time 
scales are accessible that are currently out of the reach 
of laser-based techniques. We therefore suggest that IXS, 
because of its kinematic flexibility, may provide an alter- 
native window on the attosecond phenomena targeted in 

Refs. HSU SUB II 11 

A sensible starting point for 
comparison would be the photofragmentation reaction of 
Nal or LiF. 

Not all attosecond phenomena may be probed with 
IXS. Like x-ray diffraction, it is only sensitive to phe- 
nomena which modulate the electron density, that is, 
which are capable of screening charge, such as collec- 
tive modes like plasmons and phonons. This makes IXS 
somewhat inappropriate for the study of low energy sin- 
gle particle excitations except in materials which can- 
not sustain collective electronic vibration, such as large 
gap insulators [2^. This is also true of resonance tech- 
niques which rely specifically on the coulomb interaction 
to reach final states |29, bM |31j. On the other hand, low 
energy electronic screening, which is determined by the 
plasmon through its real part, is quite accessible. 
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FIG. 3: Time frames of x( x i t) derived from Figure 1, in units 
of A -6 . The vertical scale has been clipped at 1 A~ 6 in frames 
(a-g), 0.1 A" 6 in (h), and 0.005 A~ 6 in (i). Distances are 
indicated with scale bars, (a) At t < 0, before the pertur- 
bation, the system is "placid", (b) Shortly after the initial 
event, showing a large (off scale) negative recoil at \x\ = 
surrounded by compensating positive build up. (c-g) The evo- 
lution of this disturbance at selected later times, (h) Same 
as (b) but on an expanded scale, showing distant Friedel-like 
oscillations, (i) After the disturbance has damped, scale ex- 
panded to show the experimental noise level. 



FIG. 4: Electron density disturbances for two extended 
sources, constructed via Eq. 3. (a) Induced charge around 
an oscillating dipole, such as a photoexcited chromophore in 
solution, and (b) the wake produced in water by a 9 MeV gold 
ion. The ion location and direction of motion are indicated 
by the red arrow. 



* Current address: Bldg. 725D, Brookhaven National Lab- 
oratory, Upton, NY, 11973 
[1] M. Drescher, et. al, Science, 291, 1923 (2001) 
[2] M. Drescher, et. al, Nature, 419, 803 (2002) 
[3] A. Baltuka, et. al, Nature, 421, 611 (2003) 
[4] M. Hentschel, et al, Nature, 414, 509 (2001) 
[5] R. Kienberger, et. al, Science, 297, 1144 (2002) 
[6] H. Niikura, F. Legar, R. Hasbani, M. Y. Ivanov, D. M. 

Villeneuve, and P. B. Corkum, Nature, 421, 826 (2003) 
[7] H. Niikura, F. Legar, R. Hasbani, M. Y. Ivanov, D. M. 

Villeneuve, and P. B. Corkum, Nature, 417, 917 (2002) 
[8] P. M. Paul, et. al., Science, 292, 1689 (2001) 
[9] W. Schiilke, in Handbook on Synchrotron Radiation, ed. 
D. E. Moncton and G. Brown, (North-Holland, Amster- 
dam, 1991), V. 3, p. 565 
[10] Leon van Hove, Phys. Rev., 95, 249 (1954) 
[11] K. Sturm, Z. Naturforsch., 48a, 233 (1993) 
[12] J. D. Jackson, Classical Electrodynamics, 2 ed., § 7.10 

(John Wiley & Sons, New York, 1975) 
[13] There are cases, such as near a charge density wave in- 
stability, in which the Kramers-Kronig relations are vi- 
olated at nonzero momentum transfer. See P. Martin, 
Phys. Rev., 161, 143 (1967) 
[14] There is a trick in some cases for circumventing the KK 
transform. See V. A. Madsen and G. R. Satchler, Phys. 
Rev. C, 48, 1221 (1993) 
[15] P. M. Platzman and P. A Wolff, Waves and Interactions 
in Solid State Plasmas, (Academic Press, New York and 
London, 1973) 

[16] H. Hayashi, N. Watanabe, Y. Udagawa and C.-C. Kao, 

Proc. Nat. Acad. Sci., 97, 6264 (2000) 
[17] x(k,^) is related to the refractive index by n 2 {k.,u) = 

1 - (4Tve 2 /hk 2 )x(k,u). At u = 2 eV as k -> our data 

converge to n = 1.37, close to the correct value n = 1.33. 
[18] K. D. Finkelstein, P. Abbamonte, and V. O. Kostroun, 

Proc. SPIE, 4783, 139 (2002) 
[19] J. Friedel, Phil. Mag., 43, 153 (1952) 
[20] These oscillations arise from the k dependence of x(k, w) 

in the same way that Friedel oscillations arise from the 



4 



static susceptibility, x(k). 
[21] D. Pines, Elementary Excitations in Solids, (Addison- 

Wesley, Reading, Mass., 1963) 
[22] J. Mathews and R. L. Walker, Mathematical Methods of 

Physics, (W. A. Benjamin, Inc., New York, NY, 1964), 

§9 

[23] L. Stryer and R. P. Haugland, Proc. Natl. Acad. Sci., 58, 
719 (1967) 

[24] D. R. Penn, Phys. Rev. B, 35, 482 (1987) 
[25] W. Schiilke, Phys. Lett. A, 83, 451 (1981) 
[26] J. A. Golovchcnko, el. ai, Phys. Rev. Lett., 46, 1454 
(1981) 

[27] V. M. Agranovich and V. L. Ginzburg, Spatial Disper- 
sion in Crystal Optics and the Theory of Excitons, (John 



Wiley & Sons, London, 1966) 
[28] W. A. Caliebe, J. A. Soininen, E. L. Shirley, C.-C. Kao, 

and K. Hamalainen, Phys. Rev. Lett., 84, 3907 (2000) 
[29] J. P. Hill, C.-C. Kao, W. A. L. Caliebe, M. Matsubara, A. 

Kotani, J. L. Peng, and R. L. Greene, Phys. Rev. Lett., 

80, 4967 (1998) 
[30] P. Abbamonte, C. A. Burns, E. D. Isaacs, P. M. Platz- 

man, L. L. Miller, S. W. Cheong, and M. V. Klein, Phys. 

Rev. Lett., 83, 860 (1999) 
[31] M. Z. Hasan, E. D. Isaacs, Z.-X. Shen, L. L. Miller, K. 

Tsutsui, T. Tohyama, and S. Maekawa, Science, 288, 

1811 (2000) 



This figure "ultraFigl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0311373vl 



This figure "ultraFig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0311373vl 



This figure "ultraFig3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0311373vl 



This figure "ultraFig4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0311373vl 



